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ABSTRACT 

We numerically calculate the energy and momentum transfer rates due to 
Coulomb scattering between two fluids moving with a relative velocity. The 
results are fitted by simple functions. The fitting formulae are useful to simulate 
outflows from active galactic nuclei and compact high energy sources. 



Subject headings: Relativistic thermal plasma 
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INTRODUCTION 



Relativistic jets are observed in active galactic nuclei (AGN) and Galactic black hole 
candidates. The velocity of these jets are highly relativistic with a bulk Lorentz factor above 
10. The kinetic power is almost comparable to the Eddington luminosity. The production 
and bulk acceleration of these jets are still unknown, though many ideas have been proposed 
ranging from magneto-hydrodynamical to radiative and thermal ones. There is no consensus 
on how jets are produced and accelerated. 

Although it is difficult to determine the matter c ontent of jets from observations, several 



independent arguments favor electron-positron jets flTakahara Ill994l. 



1996 


War die et al. 


1998 


Kino & Takahara 


2004 



Homan &: Wardle 111999: Hirotani et al. 



199 



i 



199 



i 



200C 



Revnolds et al 



Hirotani 



20051 : 



Croston et aLlboOsh . Electron-positron jets are most likely pro- 



duced in accretion disks around the central black holes. Because the electron mass is much 
smaller than the proton mass, the produced electron-positron pairs can be ejected more eas- 
i ly than protons. Some papers discuss the accretion disks with electron-positron outflows 



(IMisra fc Melia Ill995l : iLiang k Li Ill995l : iLi fc Liang Ill996l : lYamasaki et al. 1119991 ). If the 



accretion disks form hot pair plasma strongly coupled with photon s, the plasma may be 



thermally ac celerated like the fireball applied to gamma-ray bursts. Ilwamoto fc Takahara 



( I2OO2I . I2OO4I ) showed that a "Wien fireball" , which is optically thick to Compton scatter- 
ing but thin to absorption, results in a relativistic outfiow avoiding the difficulties of pair 
annihilation and radiation drag. 

Pairs are formed via photon-photon collisions in the accretion disk composed of normal 
plasma (electro n-proton) . Such pairs may escape from the disk by their own pressure or 



radiative force (lYamasaki et al. 1119991 ). In order to investigate pair ejection from the disk, 
one needs to treat multi-component plasma dynamics with radiation field. However, there has 
been no study of formation and ejection of pairs from the accretion disk taking into account 
radiative transfer consistently. Such studies require knowledge of friction force between the 
background plasma and pair outfiow. Although several plasma effects may be important. 
Coulomb scattering is the first one to be taken into account. The purpose of this paper is 
to obtain useful formulae of the energy and momentum transfers between two fiuids moving 
with a relative velocity via Coulomb interaction. 

The heat ing and cooling rates of thermal plasm a s are wel l known in both non-relativistic 



( ISpitzerlll956l ) and relativistic cases (IStepney 



19831 ). iDermerl (119851 ) analytically derived the 



energy exchange rate of two isotropic plasmas interacting at different temperatures. If there 
exists a relative velocity between two plasmas, analytical expressions of the energy and 
momentum exchange rates are hard to be obtained. 
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In this paper, we numerically obtain the energy and momentum exchange rates due to 
Coulomb scattering between two Maxwell-Boltzmann plas mas with a relative veloci ty. In 
another paper, we simulate pair outflows from hot plasmas (lAsano fc Takahara II2006I ) using 
the results in this paper. In §2 we describe a formulation of the reaction rates between two 
plasmas. In §3 we show the numerical results for parameter regions we are interested in. 
Fitting formulae are presented so that one can use our results in a numerical simulation of 
two fluid dynamics. 



2. TWO-FLUIDS INTERACTION 

In this section we describe a formulation of the exchange rates of energy and momentum 
due to Coulomb scattering in two fluids. Consider a collision of two particles, which belong 
to the fluids A and B, respectively. Before the collision their 4-momenta in the laboratory 
frame are 

= mAC7A(l,/3A), = mBC7B(l,/3B)- (1) 

Here, mi and /3j are the mass and velocity normalized by the light speed c of the particle in 
the fluid i, respectively. The Lorentz factor of the particle is denoted as 71 = (1 — 
and that of the relative velocity /^^ is written as 

7r = 7A7B(l-/3A-/3B), (2) 

which is Lorentz invariant. 

The velocity of the center-of-mass (CM) in the laboratory frame is 

niA'jA + mBjE 

Let this direction define the x-axis in the Cartesian coordinate {x, y, z) in the laboratory 
frame, /^q^ = I^cmx. The velocity of particles in the fluid i is decomposed as 

A = I3ixx + iSiyij + (3;^z, (4) 

Denoting the 4-momentum of particles in the fluid A after the collision as — 
"T'AC7a2(15 /3a2)) the exchange of 4-momentum is written as 

Ap^ = mAc(7A2 - 7a, 7a2/3a2 - 7a/3a) • (5) 

The mean exchange of energy for elastic scattering in the laboratory frame is 

(AE) = (Ap°)c=-(Ap°)c 



-4- 



2mAC^7A7cMsin^(a/2) /3cm(/3cm - Pax) 

(6) 



-2mA^B sin {a/2) 2 
^2 I ^2 



X [mB7A - "^a7b + 7r("^A7A - "2b7b)] 



( 1Stepneylll983l ). where a is the scattering angle in the CM-frame and 7cm = (1 ~ Pcm) 



The mean momentum exchange in the laboratory frame is described as 

(ApX) = 2mAC7A7cM(/?CM - /?Ax) sin2(a/2) 

{Apl) = -2mAC7A/?A,sin2(a/2) (7) 
(Apl) = -2mAC7A/?A.sin^(a/2) 

Hereafter we denote the distribution functions of the fluid i as where n; is the 

number density of the fluid in the laboratory frame so that J fi{Pi)d^Pi = 1. The scattering 
rate per unit volume, 

Rab = wa/a PaF^b/b Pb c — dQ, 8 

7a7b dil 



( IWeaveii Il976l : iLandau fc Lifshita Il975l ). is Lorentz-invariant. In this expression d^l = 
2tt sin ada and da/dQ is the differential cross-section. The exchange rate of 4-momentum 
per unit volume dV and per unit time dt is described as 



dP^' 1 



dtdV 1 + Sab 



j j {Ap'')RABd'pAd'pB, (9) 



Then, we now restrict the situation that particles in the fluid A have an isotropic 
Maxwell-Boltzmann (MB) distribution in the laboratory frame K and that particles in the 
fluid B have another isotropic MB distribution in another inertial frame K'. The relative 
velocity between the two frames is defined as and its Lorentz factor is Fr = (1 — 

The relativistic MB distribution of the fluid A in the laboratory frame is described as 
fA(PA)d'pA = ^^T^YU^AduAdnA, (10) 

47r6'Ai^2(l/6'A) 

where Ui = 7i/5i, 6i = Ti/mjC^ is the temperature normalized by the mass of the particles, 
K2 is the modified Bessel function of the 2nd kind, dQi = d^\d(f)\ is the solid angle of p;, 
respectively. Because of Lorentz invariance of the distribution function, ubJb = ''^b/b — 
n'^ exp{—'y'^/ 9b) / {4:Tt9bK2{1 / Ob)) , where primed values represent quantities in the frame K'. 
Since = ub/T^,, we obtain 
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where 7b = rR(7B - /^R/iBMe)- 



When the particle species in the fluids A and B are the same, we cannot determine 
if the scattered particle belongs to fluid A or B by the quantum mechanics principle. In 
this case, we treat as follows. As is well known, plasma relaxation is achieved by mainly 
small angle scattering {a ^ 1). Therefore, it may be natural to consider that a < it/ 2 in 
this case. Even if we allow large angle scattering {a > 7r/2), we ought to consider that it 
means exchange of particles between the fluids A and B. As a result, we can use the above 
formulation assuming a < 7r/2 for scatterings of the same species of particles. 

The integral over a is analytically possible (see Appendix). On the other hand, from 
the axial symmetry the sextuple integral in equation Qj is reduced to a quintuple integral. 
However, further reducti on in the order of integral may not be carried out differently from 
isotropic cases such as in iDermerl (119851 ). Therefore, a straightforward method in numerical 
integratio n is inefficient in this case. Using a Monte Carlo technique that is similar to the 
method in lRamaty fc Meszarod (jl98ll ). we numerically integrate equation (Q with N = 10^ 
collisions of two particles whose momentum distributions are proportional to equations fllOl) 
and (II ip . respectively. The differential cross-sections we used are summarized in Appendix. 
The Coulomb logarithm In A = 20 is adopted throughout this paper. 



3. NUMERICAL RESULTS 

We calculate energy and momentum transfers (ET and MT) for probable temperatures 
of leptons and protons in hot accretion disks. The relative velocity of two fluids, t/^ = 
a/F^ — 1, is in the range of 10^^ < [/r < 10^ in our calculation. Because of the property 
of the Monte Carlo method, estimated values of ET and MT may have errors especially for 
MT with [/r <^ 1. However, from our experiences, the errors are at most 10 % for ur > 0.1. 
In contrast to the uncertainty in the Coulomb logarithm, these errors may be negligible. 

We express the energy gain rate of the outflowing plasma (fluid B) as 

dE 2 / \ 

UAnBrrieC Fe{Ta,Tb,Ur). (12) 



dtdV 



in the laboratory frame (the comoving frame of the fluid A). The value has a dimension 
of [cm^/s]. We should notice that ua and ne = ^ru^ are densities in the laboratory frame. 

On the other hand, the momentum loss rate of the outflowing plasma is expressed as 

dP 

—— = nAUBrriecFpiTA, Tr, ^/r). (13) 
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The dimension of Fp is the same as that of Fe- 

Since dtdV is Lorentz invariant, these rates in the laboratory frame correlate with the 
rates in the comoving frame of the fluid B as 

dE ^ dE' dP' 

rRiirjTT; - ct/R— — (14) 



dtdV '"dt'dV "dt'dV 
dP ^ dP' Ur dE' 



(15) 



dtdV ^dt'dV c dt'dV 
respectively, where the minus signs in the second terms are due to the definition of dP/dtdV 
(momentum "loss" rate). The above equations show that negative energy gain rate {dE/dtdV < 
0) does not necessarily imply a temperature drop of the fluid B. If the fluid B is decelerated 
rapidly {dE' /dt'dV <^ cdP' /dt'dV), the energy gain rate in the laboratory frame can be 
negative. 



3.1. ENERGY TRANSFER in p-e INTERACTION 



In this subsection we show the numerical results of ET between the background protons 
(A=p) and outflowing electrons or positrons (B=e). We assume the temperature of protons 
is higher than that of the electrons (positrons). Hereafter we denote the temperature of 
protons normalized by electron mass as 0p = Tp/nicC^ = {mp/mc)Op. In Figure 1 we plot 
Fe for 6p = 10. If U-r = 0, the energy gain rate is analytically obtained as 



FE(rp,Te,o) 

2(ge + ^p)' + l 
de + 



30"t 



T - T 



2mpcK2il/e,)K2{l/ep 

6*0 + 6'p \ , ^ I Oe + 



In A X 



+ 2Ko 



(16) 



( IStepney &: Guilbertlll983l ). Our numerical results for [/r <^ 1 agree with the above analjd;- 
ical estimate. For Ur > 1.0, Fe declines with [/r. If [/r is substantially larger than 9e, we 
can neglect the thermal motion of leptons and Fe becomes independent of 6^. Therefore, Fe 
agree with each other for Up^^ 9^. 

Unfortunately, we could not fit our results by a simple analytical function of 9p, 9e, and 
[/r. In Figure 2 we plot Fe for Bp = 10° ''^ as well as fitting functions. We fit the numerical 
results by two functions of U^; exponentially damped one for lower U^, and power-law one 
for higher [/r as 



Al exp 



- m: 



for Un < f/th 
for f/R > t/th. 



(17) 
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Fig. 1. — Fe of p-e interaction for 6p = 10.0 in cgs unit. The temperatures of electrons 
(positrons), 6*6, are from 10^^ to 10° at 10°'^^ intervals in logarithmic scale. 

The energy gain rate of leptons Fe generally decreases with and for f/^ ^ 0p it becomes 
negative [see eq. ([6])]. This corresponds to the deceleration of the lepton beam moving 
through cold medium. Since ?7r ^ nip/mc in our case, /3cm ^ 1, namely the CM frame 
is almost the same as the laboratory frame. For [/r ^ Op, the Lorentz factor of electrons 
in the CM frame 7 (see APPENDIX) is approximated as ~ 7e ~ J/r ~ 7r. The equation 
indicates the energy loss of electrons in one collision (AF) ~ {mc/mp)U^meC^ sin^ ('^/2), 
while the effective cross-section (see APPENDIX) cxeff oc 7"^ oc U^'^. Therefore, Fe (roughly 
proportional to C(Tqs{AE)) becomes negative and constant for [/r ^ 9p as shown in Figure 
2. 



Even for positive Fe the fitting functions f[T7j) deviate for [/r ~ Op as shown in Figure 2. 
However, the fitting formulae may be practically correct for [/r < Gp. The fitting parameters 
^L, P, Uq, All, r, and f/th are listed in Tables 1-5. The parameter p is fixed as 2 except for 
Tables 1 and 5. We listed only for Gp = 10°■^^ 10\ 10^■^ 10^, and lO^-^. One can interpolate 



o 




0.01 0.1 1_ 10 100 



Fig. 2. — Two examples of the fitting functions for Fe of p-e interaction. Circles and 
triangles are numerical results. Filled symbols are |Fe| for negative Fe. Solid lines are the 
fitting functions we obtained (see equation f|T7|) and Table 1). 
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Fig. 3.- 



Same as Fig.l but plotted vs. ^e- 
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Fe for general values of Gp from our fitting formulae. Our fitting formulae give a satisfactory 
fit to numerical results for these ranges of 6p and [/r. 

Although the numerical results smoothly change with 6'e as shown in Figure 1, the fitting 
parameters do not necessarily change monotonically. There exist numerical uncertainties of 
~ 10 % and the fitting method is not unique. For example, if we refit the results with a slight 
change of the "connecting point " Uth, the other parameters, especially indices p, q (see next 
subsection), and r, will change substantially. Therefore, as will be shown throughout this 
paper, parameter values we obtained do not always change systematically with ^e- 

In Figure 3 we plot Fe for Op = 10 as functions of 6c for reference. The results are not 
necessarily monotonically decreasing functions of 6'e, while curves in Figure 1 decrease with 
increasing U^. 



The numerical results of Fp for Op = 10 are plotted in Figure 4. The value of Fp 
increases with increasing f/^ for f/a <^ 1. At [/r ~ 1 Fp starts to decrease because of 
smallness of the cross sections in relativistic collisions. The results are fitted in the same 
way as Fe. The fitting formulae are written as 



and the parameters are listed in Tables 6-10. As shown in Figure 5 the results are fitted by 
these functions within ~ 10 % errors for a wide range of [/r. 

Though a sufficient amount of momentum is exchanged for each collision of particles, 
the total momentum exchange between the two fluids is almost cancelled out for [/r -C ^c- 
Since Fp comes from the bulk motion of the outflowing electrons, the behaviour, Fp oc [/r 
for L^R <C 1, is a natural consequence. On the other hand, Fp decreases with [/r for [/r > 1. 
From equations (El) and ((71), the momentum loss of electrons (Ap) ~ (AF)/(c/3cm) for 
f/R ^ Op. As discussed in §3.1, (AF) oc and Pqm ?7r in this case. Then, we obtain 
Fp oc f/£^, which is consistent with the index r in the Tables. 

In Figure 6 we plot Fp for Op = 10 as functions of Og for reference. The line for 
[/r = 10~^ is not smooth probably owing to a statistical error in the Monte Carlo integral. 
Although the rate of momentum transfer is monotonically decreasing functions of 6e, the 
behavior beyond and below [/r ~ 1 is apparently different. There exists a maximum value 
of Fp at [/r ~ 1 and the Coulomb friction becomes small as 6^ increases, because of the 



3.2. MOMENTUM TRANSFER in p-e INTERACTION 




for Ur < f/th, 
for f/R > f/th, 



(18) 




0.01 0.1 1 10 100 

Fig. 4. — Fp of p-e interaction for 0p = 10.0 in cgs unit. The temperatures of electrons 
(positrons) are same as those in Fig. 1. 

monotonic behaviour of the cross section, oc 7~^/3~^ (see APPENDIX). 

3.3. ENERGY TRANSFER in e-e INTERACTION 

In this subsection wc show the resuhs for ET between the background electrons (A=bc) 
and outflowing electrons (B= — ) or positrons (B= +). We assume that the temperature 
of the outflowing leptons is lower than the background electrons. In the opposite cases ET 
and MT are obtained from the Lorentz transformation of our results. The values for e~-e~ 
and e~-e+ interactions are almost the same. Therefore, we only show the results for e~-e~ 
hereafter. 

Figure 7 shows ET between two electron fluids for = 1- For 6_ < 0.1 the functional 
shapes of ET are almost invariant, because we can neglect the bulk motion for C/r <g ^be- 
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I I I I I I I I I I I 

0.01 0.1 10 100 

Fig. 5. — Two examples of the fitting functions for Fp of p-e interaction. Circles and 
triangles are our numerical results. Solid lines are the fitting functions we obtained (see 
equation (fT8i) and Table 5). 
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Fig. 6. — Same as Fig.4 but plotted vs. 9c- Dashed lines are from [/r =10 ^ to 10 and 
solid lines are from U^. = 10^ to 10^ at 10°'^ intervals in logarithmic scale. 




Fig. 7. — -Fe of e^-e^ interaction for 6'be = 1 in cgs unit. The temperatures of the outflowing 
electrons, 6'_, are from 10~^ to 10*^ at 10°'^^ intervals in logarithmic scale. Dashed lines are 
|Fe| for negative F^. 

From the figure we can see ET becomes zero at f/^ ~ a few, although it is difficult to 
determine the value of f/a, at which ET becomes zero, within our computational precision. 
As 6'_ increases, |Fe| also increases in cases of negative ET ([/r > 1), while |Fe| decreases 
for positive ET (f/^ < 1). 

When ^ ^^be, the thermal motion can be neglected. In this case, from equation ([6]), 
the energy loss of outflowing electrons in one collision {AE) ~ U^m^c^ sin^ (c^/S)- Since the 
mass of the colliding particle is the same, the Lorentz factor of electrons in the CM frame 
7 ~ 7cM U^'^ . The effective cross section crgfr oc 7"^ oc U^^ so that Fe becomes constant 
for [/r ^ 9y,c as shown in Figure 7. 

We fit our results by functions. 



Fe = ^Lexp [-[U^/U^f], 



(19) 



0.01 



0.1 



e 



Fig. 8. — Fe of e -e interaction for L^r = in cgs unit. The temperatures of the background 
electrons, ^bej are from 10~^ to 10°-^ at 10°-^^ intervals in logarithmic scale. 



for C/r < 1 and 

Fe = -^H exp [- {Un/Unf], (20) 
for Uji>l. Fitting results are listed in Tables 11 and 12. 

We tabulate the fitting parameters only for 6}^^ = 0.1 and 1. However, as is shown in 
Figure 7, Fe is almost constant for Ur <^ 1. Thus, it is valuable to estimate Fe for = 
with a wide range of ^be- Figure 8 shows the results, which are fitted by a function 

FE = Aexp[-(^_/^or], (21) 

(see Table 13). Although the 6'be-dependece is hardly expressed by a simple formula, we can 
interpolate Fe from Table 13. 
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Fig. 9. — Fp of e -e interaction in cgs unit. 



3.4. MOMENTUM TRANSFER in e-e INTERACTION 

As is the case with energy transfer, the momentum transfers for e^-e^ and e^-e+ interac- 
tions are almost the same. The numerical results of Fp are plotted in Figure 9. For = 10^^ 
the results are almost independent of the temperature of the outflowing electrons Even 
for ^be = 10~^, Fp decreases only by a factor of less than 2 with increasing ^_ = 10~^ to 10~^. 
Figure 9 clearly shows the reduction of Fp due to increasing ^be- As we have seen in the 
former cases, the growth of the average 7 by rising ^be results in the reduction of the cross 
section ((Jeff oc 7"^/?"^ for relativistic case). For [/r ^ 1, Fp is almost constant differently 
from the case of p-e interaction. This is because /9cm ~ 1 in e-e interaction and a resultant 
relation (Ap) ~ (AF)/c. Therefore, our results show Fp ~ — Fe oc t/^ for t/R ^ 6'be- 

In contrast to p-e interaction, the cooling and heating properties of the two fluids in this 
case are the same. Therefore, if the flow velocity is non- relativistic, the two temperatures 
may be the same. Thus we calculate Fp mainly with ^_ — O^e, as shown in Figure 10. Even 




0.01 0.1 Itt 10 100 

Fig. 10. — Fp of e^-e" interaction for ^bc = in cgs unit. The temperatures are from 10~^ 
to lO"'^ at lO"'^^ intervals in logarithmic scale. 

if the two temperatures are different, the correction is not so large as shown in Figure 9. We 
fit our results by functions; 




fort/R>t/th, 



(22) 



and results are summarized in Table 14. The symbol oo for Un in this table means that Fp 
is constant as Fp = An for C/r > Uth- 




Fig. 11. — Fp of e -e"*" interaction due to scattering (dashed) and pair-annihilation (sohd) 
for ^be = 10° in cgs unit. 

3.5. MOMENTUM TRANSFER in PAIR-ANNIHILATION 

We shortly comment on the effect of pair-annihilation process here. The outflowing 
positrons may interact with the background electrons and turn into gamma-rays. This 
process decreases the momentum of the outflowing fluid, and it is not trivial whether we 
can neglect this momentum loss process or not. In Figure 11 we plot Fp of e^-e"*" due to 
pair-annihilation. As ^_ increases for ^^e — 10°, the momentum loss due to pair-annihilation 
increases, while that due to pair-scattering is almost constant. Even for 9- — 10°, the 
momentum loss due to pair- annihilation is smaller than that due to pair-scattering by one 
order of magnitude. 

For smaller values of ^bc < 10°, the contribution of pair-annihilation to Fj> becomes 
much more negligible. Therefore, we can neglect the effect of pair-annihilation on the flow 
within precision we require. 
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4. SUMMARY 

Motivated by electron-positron outflows from AGNs, we numerically calculate the energy 
and momentum transfer rates due to Coulomb scattering between two fluids with a relative 
velocity for plausible parameters in AGN models. Although several plasma effects, such as 
two-stream instability etc., may enhance the energy and momentum transfers, the effects 
of Coulomb scattering are the first to be taken into account. Our tables obtained from 
the numerical results are useful to simulate pair outflows from hot plasmas, or evaluate the 
interaction between AGN jets and the ambient medium. Especially, the momentum transfer 
ra te is indispensable f o r suc h simulations. Using the results in this paper, a simulation 



m 



Asano fc Takahara I (120061 ) shows that the frictional force due to Coulomb scattering is 
comparable to radiative force for plausible parameter sets in AGN jet models. Therefore, we 
cannot neglect the effects of Coulomb scattering in such simulations. 

We thank the anonymous referee for his useful comments. This work is partially sup- 
ported by Scientific Research Grants (F.T. 14079205 and 16540215) from the Ministry of 
Education, Culture, Science and Technology of Japan. 
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Table 1: Fitting parameters of for p— e interaction, where Gp = 10° ''^. See eq. (1171) . 
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1.443 


1.3 








IQO.O 


1.7e-15 


1.511 


1.3 









Table 2: Same as Table 1 but for Op = 10 



00 


Al 


Uo 


An 


r 


f/th 


10-2 


1.4e- 


12 


0.179 


5.6e- 


•15 


-1.2 


0.37 


10-1.75 


7.3e- 


13 


0.224 


5.6e- 


•15 


-1.2 


0.44 


10-1.5 


3.7e- 


13 


0.284 


5.8e- 


•15 


-1.2 


0.52 


10-1.25 


1.7e- 


13 


0.374 


5.9e- 


•15 


-1.2 


0.63 


10-1 


7.7e- 


14 


0.515 


6.2e- 


•15 


-1.2 


0.76 


10-0.75 


3.5e- 


14 


0.698 


7.7e- 


•15 


-1.4 


0.69 


10-0.5 


1.5e- 


14 


0.909 


5.2e- 


•15 


-1.4 


0.76 


10-0.25 


7.0e- 


15 


1.161 


3.6e- 


•15 


-1.4 


1.01 


10° 


3.2e- 


15 


1.544 


2.3e- 


•15 


-1.4 


1.31 


100.25 


1.6e- 


15 


1.359 


9.8e- 


•16 


-1.4 


1.16 


100.5 


7.1e- 


16 


1.209 
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Table 3: Same as Table 1 but for Op = 10^-^ 



0e 








Ah 


r 




Uth 


10- 


-2 


i.oe-i.^ 




1 1 A 

i.ye-i4 


-i 


.z 


U.40 


10- 


-1.75 


1.3e-12 


0.274 


2.0e-14 


-1 


.2 


0.50 


10- 


-1.5 


7.7^13 


0.303 


2.0^14 


-1 


.2 


0.51 


10- 


-1.25 


4.1e-13 


0.425 


2.3e-14 


-1 


.2 


0.66 


10- 


-1 


2.1e-13 


0.545 


2.2e-14 


-1 


.2 


0.75 


10- 


-0.75 


9.9e-14 


0.778 


2.4e-14 


-1 


.3 


0.87 


10- 


-0.5 


4.7e-14 


0.952 


1.7e-14 


-1 


.3 


0.79 


10- 


-0.25 


2.2e-14 


1.230 


1.3e-14 


-1, 


.3 


1.04 




l.le-14 


1.479 


6.9e-15 


-1 


.2 


0.96 


100.25 


5.6e-15 


1.766 


4.6e-15 


-1 


.3 


1.19 


100.5 


2.9e-15 


1.783 


2.5e-15 


-1 


.3 


1.42 



Table 4: Same as Table 1 but for Gp = 10^ 



0e 






Uo 


Au 


r 


C/th 


10- 


■2 


1.6e-12 


0.383 


6.1e-14 


-1 


0.65 


10- 


■1.75 


1.4e-12 


0.408 


6.2e-14 


-1.1 


0.66 


10- 


■1.5 


l.le-12 


0.446 


6.4e-14 


-1.1 


0.69 


10- 


■1.25 


7.3e-13 


0.524 


6.9e-14 


-1.2 


0.74 


10- 


■1 


4.4C-13 


0.644 


7.3e-14 


-1.2 


0.80 


10- 


■0.75 


2.5C-13 


0.825 


6.0e-14 


-1.2 


0.96 


10- 


■0.5 


1.3e-13 


0.933 


5.1e-14 


-1.2 


0.75 


10- 


0.25 


6.4e-14 


1.245 


3.5e-14 


-1.1 


0.95 


10° 




3.3e-14 


1.441 


2.2e-14 


-1.1 


1.10 


100.25 


1.7e-14 


1.630 


1.3e-14 


-1.1 


1.23 


10" 




9.1e-15 


1.816 


7. Be- 15 


-1.1 


1.37 
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Table 5: Same as Table 1 but for Op = lO^-^ 



Be 






Uo 


P 


An 


r 




Uth 


10-2 


l.le-12 





.673 


2.0 


2.2e-13 


-1, 


,1 


0.79 


10-1.75 


l.Oe-12 





.704 


2.0 


2.1e-13 


-1. 


,1 


0.83 


10-1.5 


9.4e-13 





.732 


2.0 


2.1e-13 


-1. 


,1 


0.83 


10-1.25 


8.0e-13 





.782 


2.0 


2.1e-13 


-1. 


,1 


0.84 


10-1 


6.2e-13 





.867 


2.0 


2.0e-13 


-1. 


,1 


0.86 


10-0.75 


4.3e-13 





.992 


2.0 


1.7e-13 


-1. 


,1 


0.87 


10-0.5 


2.7e-13 


1 


.164 


2.0 


1.3e-13 


-1. 


,1 


1.05 


10-0.25 


1.5e-13 


1 


.540 


1.6 


9.5e-14 


-1, 


,1 


1.36 


10° 


8.5e-14 


1 


.757 


1.6 


6.5e-14 


-1. 


,1 


1.45 


100.25 


4.7e-14 


1, 


.924 


1.6 


4.1e-14 


-1. 


,1 


1.73 


100.5 


2.6e-14 


1, 


.994 


1.6 


2.4e-14 


-1, 


,1 


1.68 



parameters of Fp 


for p— e interaction 


, where B 


p 


100.75_ , 


ee 


Al 


Uo 


P 


q 


An 


r 


f/th 


10-2 


l.le-12 


0.214 


1.5 


1.0 


8.7e-13 


-1.1 


0.48 


10-1.75 


4.6e-13 


0.255 


1.5 


1.1 


7.8e-13 


-1.1 


0.56 


10-1.5 


2.2e-13 


0.322 


1.5 


1.1 


7.8e-13 


-1.1 


0.65 


10-1.25 


l.le-13 


0.502 


1.5 


1.0 


7.8e-13 


-1.1 


0.90 


10-1 


5.0e-14 


0.631 


1.5 


1.0 


5.6e-13 


-1.0 


1.16 


10-0.75 


1.7e-14 


0.700 


1.3 


1.1 


5.6e-13 


-1.0 


1.02 


10-0.5 


5.9e-15 


0.713 


1.0 


1.2 


4.1e-13 


-1.0 


2.02 


10-0.25 


4.0e-15 


0.572 


0.8 


1.2 


3.1e-13 


-1.0 


2.29 


lOO.o 


2.3e-15 


0.798 


0.8 


1.1 


2.2e-13 


-1.0 


2.70 



See eq. ([T8]). 
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Table 7: Same as Table 6 but for Bp = 10 









TT 


P 




A 


r 


TT 


10- 


-2 


8.6e-13 


0.237 


1 


.5 


1.0 


l.le-12 


-1.2 


0.44 


10- 


-1.75 


4.0e-13 


0.257 


1 


.5 


1.1 


8.1e-13 


-1.1 


0.53 


10- 


-1.5 


1.9e-13 


0.341 


1, 


.5 


1.1 


7.9e-13 


-1.1 


0.66 


10- 


-1.25 


l.le-13 


0.488 


1, 


,5 


1.0 


7.5e-13 


-1.1 


0.85 


10- 


-1 


3.1e-14 


0.512 


1, 


,5 


1.2 


6.8e-13 


-1.1 


0.74 


10- 


-0.75 


1.8e-14 


0.914 


1, 


.5 


1.0 


5.6e-13 


-1.0 


0.88 


10- 


-0.5 


9.5e-15 


1.074 


1, 


.5 


1.0 


4.1e-13 


-1.0 


1.52 


10- 


-0.25 


4.7e-15 


1.195 


1 


.2 


1.0 


3.0e-13 


-1.0 


1.86 


1000 


1.8e-15 


1.090 


1 


.0 


1.1 


1.7e-13 


-0.9 


2.75 


100.25 


1.2e-15 


1.418 


1 


.0 


1.0 


1.2e-13 


-0.9 


3.23 


100.5 


9.9e-16 


1.617 


1 


.0 


0.9 


6.9e-14 


-0.9 


3.44 



Table 8: Same as Table 6 but for Bp = 10^'^ 







Uo 


P 


9 




r 




10-2 


3.6e-13 


0.329 


1.5 


1.0 


7.7e-13 


-1.1 


0.69 


10-1.75 


2.5e-13 


0.359 


1.5 


1.0 


7.6e-13 


-1.1 


0.70 


10-1.5 


1.4e-13 


0.359 


1.5 


1.1 


7.5e-13 


-1.1 


0.64 


10-1.25 


8.1e-14 


0.538 


1.5 


1.0 


7.2C-13 


-1.1 


0.92 


10-1 


4.0C-14 


0.676 


1.5 


1.0 


6.6C-13 


-1.1 


1.00 


10-0.75 


1.9c- 14 


0.866 


1.5 


1.0 


5.3C-13 


-1.0 


1.15 


10-0.5 


l.lc-14 


1.271 


1.8 


0.9 


4.0C-13 


-1.0 


1.32 


10-0.25 


4.6e-15 


1.187 


1.2 


1.0 


2.9e-13 


-1.0 


1.82 


lOO.o 


2.0e-15 


1.600 


1.2 


1.0 


2.1e-13 


-1.0 


3.24 


100.25 


1.2e-15 


1.560 


1.2 


1.0 


l.le-13 


-0.9 


2.29 


100..-. 


5.5e-16 


1.191 


1.0 


1.1 


6.80-11 


-0.9 


2.39 
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Table 9: Same as Table 6 but for Bp = 10^ 





Al 


T T 

Uo 


P 


Q 


Au 


r 


T T 

Uth 


10-2 


l.Oe-13 


0.478 


1.5 


1.0 


7.2e-13 


-1.1 


0.81 


10-1-75 


7.7e-14 


0.527 


1.5 


1.0 


7.1e-13 


-1.1 


0.84 


10-1.5 


5.4^14 


0.463 


1.5 


1.1 


6.7^13 


-1.1 


0.61 


10-1.25 


5. 7c- 14 


0.697 


1.5 


0.9 


5.3e-13 


-1.0 


0.75 


10-1 


2.0e-14 


0.677 


1.5 


1.1 


5.7e-13 


-1.0 


0.85 


10-0.75 


1.5e-14 


0.876 


1.3 


1.0 


4.9e-13 


-1.0 


1.23 


10-0.5 


8.0C-15 


1.047 


1.3 


1.0 


3.6C-13 


-1.0 


1.56 


10-0.25 


4.0e-15 


1.229 


1.3 


1.0 


2.6e-13 


-1.0 


1.74 


lOO.o 


2.3e-15 


1.255 


1.0 


1.0 


1.9e-13 


-1.0 


2.51 


100.25 


8.7e-16 


1.139 


1.0 


1.1 


9.5e-14 


-0.9 


2.59 


100.5 


6.7e-16 


1.358 


1.0 


1.0 


6.3e-14 


-0.9 


2.63 



Table 10: Same as Tabic 6 but for Qp = lO^-^ 



0e 




A^ 


Uo 


P 


Q 


An 


r 




10' 


-2 


1.8e-14 


0.783 


1.5 


1.0 


4.6e-13 


-1.0 


0.97 


10- 


-1.75 


1.9e-14 


1.008 


2.0 


0.9 


4.5e-13 


-1.0 


0.99 


10- 


-1.5 


1.3e-14 


0.925 


2.0 


1.0 


4.3e-13 


-1.0 


0.94 


10" 


-1.25 


1.4c- 14 


1.102 


2.5 


0.9 


4.0e-13 


-1.0 


0.99 


10- 


-1 


9.1C-15 


1.042 


1.5 


1.0 


4.0C-13 


-1.0 


1.29 


10- 


-0.75 


6.5C-15 


1.178 


1.5 


1.0 


3.6C-13 


-1.0 


1.45 


10- 


-0.5 


4.5e-15 


1.218 


1.5 


1.0 


2.8C-13 


-1.0 


1.51 


10- 


-0.25 


2.8e-15 


1.298 


1.1 


1.0 


2.2e-13 


-1.0 


2.57 


100.0 


1.6e-15 


1.436 


1.2 


1.0 


1.3e-13 


-0.9 


2.11 


100.25 


7.3e-16 


1.124 


1.0 


1.1 


8.5e-14 


-0.9 


2.32 


10" 




1.7e-16 


0.916 


0.8 


1.1 


5.1e-ll 


-0.9 


2.89 
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Table 11: Fitting parameters of for e— e interaction, where ^be = 10 ^. See eqs. ( IT9l) and 

(EOD. 

9- Ai^ Ui^ An Uu P 

10- 2 1.2e-12 0.315 5.2e-13 0.5 -3 
10^^-^^ 9.7e-13 0.313 5.2e-13 0.5 -3 
10-^-5 6.8e-13 0.321 5.2e-13 0.5 -3 
10-^-25 3,4e.i3 0.265 5.2e-13 0.5 -3 



Table 12: Same as Table 11 but for ^bc = 10' 



e_ 






Ul 


An 


Un 


P 


10- 


-2 


6.5e-13 


1.117 


2.2e-13 


3.469 


-2 


10- 


-1.75 


6.3e-13 


1.124 


2.2e-13 


3.444 


-2 


10- 


-1.5 


6.1e-13 


1.123 


2.2e-13 


3.490 


-2 


10- 


-1.25 


5.7e-13 


1.112 


2.2e-13 


3.218 


-1.5 


10- 


-1 


5.0e-13 


1.104 


2.2e-13 


2.866 


-1.5 


10- 


-0.75 


4.0e-13 


1.076 


2.2e-13 


2.881 


-1.5 


10- 


-0.5 


2.7e-13 


0.915 


2.2e-13 


2.408 


-1.5 


10- 


-0.25 


1.3e-13 


0.679 


2.2e-13 


1.805 


-1.5 



Table 13: Fitting parameters of Fe for e— e interaction for U-r = 0. See eq. ( l2li) . 



^bc 


A 


Oo 


P 


10-1 


1.2e-12 


0.049 


1.5 


10-0.75 


l.Oe-12 


0.080 


1.4 


10-0.5 


8.5e-13 


0.13 


1.3 


10-0.25 


7.2e-13 


0.22 


1.3 


lOO.o 


6.5e-13 


0.35 


1.2 


100.25 


6.1e-13 


0.52 


1.0 


100.5 


5.9e-13 


0.7 


1.0 
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Table 14: Fitting parameters of Fp for e— e interaction for 6^ = ^bc- See eq. (1211) . 









P 


Q 


An 


Uu 


r 




10-2 


l.le-12 


0.273 


2.0 


1.0 


5.9e-13 


1.111 


-0.85 


0.36 


10-1.75 


5.3e-13 


0.428 


2.0 


0.9 


5.9e-13 


1.093 


-0.9 


0.55 


10-1.5 


1.9e-13 


0.511 


2.0 


1.0 


5.7e-13 


1.023 


-1.0 


1.16 


10-1.25 


8.9e-14 


0.676 


1.5 


1.0 


5.5e-13 


1.059 


-1.0 


0.89 


10-1 


5.0e-14 


1.081 


1.5 


0.9 


5.2e-13 


0.888 


-1.0 


1.38 


10-0.75 


1.3e-14 


0.650 


0.75 


1.2 


4.7e-13 


0.762 


-1.0 


2.33 


10-0.5 


8.8e-15 


0.931 


0.7 


1.1 


4.0e-13 


0.525 


-1.0 


3.13 


10-0.25 


3.0e-15 


0.638 


0.7 


1.3 


3.1e-13 


oo 




2.22 


lOO.o 


2.7e-15 


1.262 


0.7 


1.1 


2.2e-13 


oc 




3.27 


100.25 


1.5e-15 


1.357 


0.7 


1.1 


1.4e-13 


oo 




2.59 


100.5 


9.8e-16 


2.369 


1.0 


1.0 


8.8e-14 


oo 




2.37 
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A. DIFFERENTIAL CROSS-SECTIONS FOR COULOMB SCATTERING 

A.l. Electron-proton scattering: Rutherford cross-section 

Differential cross section is given by 

da Sot A , '2'jme \ 1 , . 

^TTTIP, ajRutherford - J-^T ^ H 77^) l^-LJ 

ail 327r7^p^ \ '^p / sm (a/2) 



( IStepney fc Guilbertlll983l ). where /3 and 7 is the velocity of electron and its Lorentz factor 
in the CM frame and ctt is the Thomson cross section. From equations (E]) and (JTj), we need 
to integrate 2sin^(a/2)^ over the scattering angles to obtain the energy and momentum 
transfer. Therefore, we define the effective cross-section as 

acs = / (1 — cosa)— 27rsina(ia. (A2) 

Here amin is the lower limit of scattering angle. Then we obtain 

3(Tt / 27me\ 

0"efr,Ruthcrford — r, o nA I "I miV, (^Adj 



272/54 \^ rUp 

where 



\nA = \n{—-^-—-}, (A4) 



sin(amin/2) 

is the Coulomb logarithm and usually evaluated about 10 ~ 20. 



A. 2. Electron-electron scattering: M0ller cross-section 



The differential cross-section of e -e scattering (jjauch fc RohrlichlllQSOl ) is given by 



— (/?, a)M0ller 



3(Tt 



^7172(72 - ly 



(27^ 



sin"^ a 



27^-7^-1^(71 



sin^ a 



(A5) 



Since the particles are not identical, a is from Omin to tt/2. In the same way as Rutherford 
scattering, we obtain the effective cross-section with A ;» 1 as 



f'"eff,M0ller 



3(Jt 

1676/54 



(272 - 1)2 [2 In A + 1] - 4 ( 27^ - 7^ 



In 2 + 



• (A6) 
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A. 3. Electron-positron scattering: Bhabha cross-section 



The differential cross-section of e -e^ scattering (jjauch fc RohrlichlllQSOl ) is given by 



da 



{13, a] 



Bhabha 



ScTx J 1 

1287172 ■ \[7/3sin(a/2)]4 
1 



+ (27/3 cos I) + 2(7/5^ (1 + cos^ I) 
3 + 2(27/3cos|)' + i(27/3cos|) 



[72/?sin(a/2)]2 
+1 [3 + 4(7/3)2 + (7/?r(l + cos2a)]|. 

The effective cross-section is given by 
3(Tt 



Ccff,Bhabha 



16 



2(272-1)2^^^ 472 + 3 67^ + 472-1 4/32-2/54 



76/34 



76/32 



276 



37^ 



(A7) 



